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National Phase of PCT/FR99/02076 
International Filing Date: September 1, 1999 

Inventors: Nathalie Laurent-Chatenet, Patrick Lechat and Henri Sanson 

Title: Method for Estimating Movement Between Two Images 

Priority: French Application No. 98 11227; Filed September 4, 1998 

Attorney Docket: 136.157 Customer No. 023907 

PRELIMINARY AMENDMENT 

DO/EO/US 

Assistant Commissioner for Patents 
Washington DC 20231 

Sir: 

This Preluninary Amendment is directed to a new U.S. application as identified above. 
Please enter this preliminary amendment prior to calculating the fees. 
Please amend the application as follows: 
IN THE SPECIFICATION 

Page 1, after the title insert the heading - BACKGROUND OF THE 
INVENTION -; 

Page 3, line 14, insert the heading -- OBJECTS AND SUMMARY OF THE 
INVENTION -; 

Page 4, line 31, insert the heading - BRIEF DESCRIPTION OF THE DRAWINGS --; 
Page 5, line 8, insert the heading - DESCRIPTION OF THE PREFERRED 
EMBODIMENTS --. 



Preliminary Amendment - Laurent-Chatenet et al. 
Attorney Docket 136.157 
Page 2 

IN THE CLAIMS 

Claims 4 and 5, line 1, cancel "one of claims 1 to 3" and substitute ~ claim 1 -; 
Claim 6, line 1, cancel "one of claims 1 to 5" and substitute - claim 1 ~; 
Claims 9, 10 and 11, line 1, cancel "one of the preceding claims" and 
substitute — claim 1 ~; 

Claims 12 and 13, line 1, cancel "one of claims 1 to 11" and substitute ~ claim 1 --; 
Claim 14, line 1, cancel "one of the preceding claims" and substitute - claim 1 -; 
Claim 16, line 1, cancel "or 15"; 

Claim 18, line 1, cancel "one of claims 12, 13 or 17" and substitute -- claim 1 --; 
Claim 19, line 2, cancel "one of claims 1 to 18" and substitute ~ claim 1 ~; 
Claim 21, line 2, cancel "one of claims 1 to 18" and substitute — claim 1 — ; 
Claim 23, line 1, cancel "one of the preceding claims" and substimte ~ claim 1 ~. 

IN THE ABSTRACT 

Please add page 21 containing an Abstract of the Disclosure. 

REMARKS 

This application has been amended to msert headings in the specification, to eliminate 
the multiple dependencies in the claims, and to add an Abstract of the Disclosure. Entry of the 
amendments and early consideration and allowance are respectfully requested. 
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ABSTRACT OF THE DISCLOSURE 
A method for estimating motion between two digital images Ii and I2, with brightness 
Yi and Y2 for generating, for each point of coordinates x, y of image h a motion vector 
(f(x,y)=(dx,dy) so as to form an image I2 from image Ii with brightness Y2(x,y)=Yi(x-dx,y-dy) 
which is an approximation of raiage I2. The method consists in using a model of finite 
elements and a differential method for determining the motion field of image 12. The mesh 
associated with the model of finite elements is refined as the motion estimation is being carried 
out. 
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METHOD FOR ESTIMATING MOVEMENT BETWEEN TWO IMAGES 

The present invention concerns a metinod for estimating movement between two 
numerical Images. 

5 

The movement between two successive images, li and b, is generally defined In 
the form of a movement field associated with either of the Images , li and b and 
constituted by movement vectors, each relating to one point of the image 
concerned. The movement vector is a two-dimensional vector representative of 
10 the difference of position between the pixel of the image h and the associated 
pixel of the Image b relating to the same physical point of the filmed scene. 

An evaluation of movement is useful In fields for processing the Image requiring a 
knowledge of movements or disparities between two images. By way of 
15 examples, It is possible to cite the following spheres of applications : 

- the compression of Images : the evaluation method Is used to limit the amount of 
data to code an Image, the images being defined In relation to one another ; 

- the compression of data In spaces of dimension greater than 2 ; 

- video coding : the movement field defined from already coded Images is then 
20 used to predict the next image ; 

- medical imagery: the method for estimating the movement between two Images 
Is used to conduct an analysis of the movement of the heart for example ; 

- remote monitoring : the method can be used to monitor road traffic ; 

- three-dimensional reconstruction from multi-view images: the method Is used to 
25 estimate disparities between various views. 

So as to obtain this movement field, a method is known on how to break down the 
image into finished elements. These finished elements, which may for example be 
triangles or quadrangles, are determined by a meshwork whose nodes 
30 correspond to the tops of the finished elements. A movement vector is calculated 
for each node of the meshwork. Then, via the bias of an Interpolation function, it is 
possible to deduce from this a movement vector for each point of the Image In 
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question. The movement field is thus determined by a model of finished elements 
defining the meshwork used to partition the image into finished elements and the 
interpolation function making it possible to calculate the movement vector at any 
point of the image. 

5 

The meshwork used can be regular or irregular and needs to selected as 
sufficiently dense so as to model as best as possible the movement between the 
two images without however requiring an excessive quantity of calculation or data 
to be transmitted. This choice is made once only at the start of the method and 
10 this meshwork generally remains the same throughout the estimate. 

The calculation of the movement vectors of the nodes of the meshwork can be 
carried out according to various methods. First of all there are putting into 
correspondence or "matching' methods consisting of testing a discrete set of 

15 possible values of movement vectors for each node of the meshwork and of 
retaining the best vectors according to a given criterion. A second method known 
as a transformed method consists of using the properties of the Fourier transform 
and its extensions so as to convert the movement into a phase jump in the 
transformed space. Finally, there is a third method known as a differential method 

20 for determining the movement vectors by optimising a mathematical criterion (for 
example a quadratic error between the image and its aforesaid value with the 
movement field). This method is most frequently used for estimating movement 
with modelisation by finished elements. A conventional differential method for 
optimising movement vectors is the Gauss-Newton method. The present 

25 application concerns more particularly the movement estimation method family 
using a model of finished elements and a differential method to determine the 
movement field. 

Although widely used, this type of method does however have several drawbacks. 
30 The meshwork selected at the start of the method may prove to be inappropriate 
with respect to the semantic contents of the image, certain zones of the image 
requiring a denser meshing and others a more aerated meshing. In addition, 



under the effect of the field of the movement vectors of the nodes of the 
meshwork, the initial meshing on the start image, for example b, is transformed 
into a new meshing on the other image, for example li. Then pathological 
situations may then occur at the level of the new meshing, such as : 

- reversals of finished elements : finished elements reverse and cover others, thus 
destroying the property of partitioning of the range of the image any meshwork 
needs to verify, 

- overflowing of the peripheral nodes of the moved meshwork after applying 
movement vectors beyond the range of the image h : certain pixels of the image I2 
can be associated with pixels of the image h situated outside the range of the 
image h. This is not strictly a problem, but it may be advantageous to force the 
peripheral nodes of the moved meshwork to remain inside the range of the image 
I1. 

One object of the invention concerns a method for estimating movement in which 
the meshwork is optimised during estimation so as to obtain at the end of the 
method a meshwork adapted to the semantic contents of the images. To this 
effect, the finished elements are refined during the movement estimate. 

Another object of the invention is to improve the effectiveness of the Gauss 
Newton method so as to optimise the movement vectors of the nodes of the 
meshwork. To this effect, this optimisation is carried out on several resolution 
levels of the images. 

Finally, another aim of the invention is to provide a method for estimating 
movement so as to avoid aforesaid pathological situations. To this effect, the 
invention provides adding during the movement vectors optimisation step 
constraints so as to avoid these situations. 

Also, the invention concerns a method for estimating the movement between two 
numerical images I1 and b with luminance Yi and Y2 and intended to generate for 
each point of coordinates x,y of the image b a movement vector d(x,y)=(dx,dy) so 
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as to form an image I2 from the image h with luminance Y2(x,y)=Yi(x-cix,y-dy) 
which is an approximation of the image b, characterised in that it comprises the 
following steps : 

(a) defining an initial model of finished elements comprising a meshing whose 
5 nodes are points of the image I2, a movement vector at each node of said 

meshing, and an interpolation formula for calculating the value of the movement 
vector of each point of the image b from the values of the movement vectors of 
the nodes of the mesh to which it belongs, 

(b) optimising the value of the movement vectors of the model according to a 
10 differential method, 

(c) calculating a variation E between the image b and the image I2 for each 
finished element or mesh, 

(d) carrying out a finer meshing on a discrete fraction of the set of finished 
elements determined according to a criterion relating to the variations E and 

15 allocating a movement vector to each new meshing node, 

(e) repeating the steps (b), (c) and (d) on the model of finished elements 
obtained at the end of the preceding step (d) until a stoppage criterion is satisfied. 

According to an improved embodiment, for each numerical image I1 and b in 
20 addition a set of R images \] is defined with a level of resolution r and luminance 
Y/ with r taking the values (0,...,R-1) and i the values 1 and 2, the images 1° and 
1° corresponding to the numerical images I1 and b, the steps (b) to (e) being 
conducted for each level of resolution r from the level r=R-1 to the level r=0. 

25 Finally, according to a preferred embodiment, constraints are added to the 
movement of the finished elements at the time of optimising the movement 
vectors so as to avoid the reversal of the finished elements. According to another 
embodiment, it is also possible to introduce constraints so as to avoid the flowing 
over of the meshing obtained after applying the movement vectors beyond the 

30 sphere of the image I1 . 

Other characteristics and advantages of the invention shall appear on a reading of 
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the following detailed description with reference to the accompanying drawings on 
which : 

- figure 1 represents a diagram of a first embodiment of the movement estimation 
method of the invention ; 

5 - figure 2 shows the step (d) of the method of the invention, and 

- figure 3 shows a diagram of an improved embodiment of the movement 
estimation method of the invention. 

Reference is first made to two numerical images h and b with respective 
10 luminance Yi and Y2. The method of the invention consists of generating for each 
point P of coordinates (x, y) in the image I2 a movement vector d(x,y)=(dx,dy). 
This vector is defined as being the vector able to construct from the image U an 
image I2 with luminance (x,y)=Yi(x-dx,y-dy) which is an approximation of I2. 
The movements are thus defined from the image U towards I2. 

15 

The sought-after movement field is defined by a model of finished elements. In 
the remainder of the description, the finished elements are regarded to be 
triangles without there being any limitation of the extent of the present application 
to this form of finished elements. As a result, the model of finished elements 
20 comprises a triangular meshing, movement vectors defined in the meshing nodes, 
and an interpolation formula for calculating the movement vector of the points 
inside the triangles. 

The interpolation formula used to calculate the movement field at any point of the 
25 range of the image I2 is the following : 

If the point P of coordinates {x,y) is considered in the image I2 belonging to the 
triangle e with vertices Pj, Pj and Pk with respective coordinates (Xj.yj), (Xj,yj) and 
(Xk,yk), its movement vector is equal to 

30 d{x,y)^J^Xix,y).d{xM 

where represents a basic function associated with the triangle e 



In the case of an affine interpolation, the '¥^^{x,y) represent the barycentric 
coordinates of the point P in the triangle e with vertices Pj, Pj, Pk. These functions 
are defined by the following equation : 



The affine functions ^f(x,y) and %^(x,y) deduced from the function ^i^(x,y) by 
circularly permuting the indices i,j,k. It is also possible to use more evolved 
models of finished elements, the functions v|/ then being able to be extended to 
polynomials with degree n=2, but the interpolation formula of the movement 
vectors then introduces first, second derivatives, etc. A miscellaneous choice of 
models of finished elements is shown in the work "Handbook of Numerical 
Analysis" by P.G. Ciarlet and J.L Lions, Volume 2, pages 59-99, published by 
North Holland. 

According to the invention, as the movement is gradually estimated, the value of 
the movement vectors of the meshing nodes known as nodal vectors is optimised 
and the meshing is locally densified when this is necessary. Advantageously, this 
optimisation shall be carried out on several resolution levels starting with a low 
resolution level. 

According to a first embodiment shown on figure 1, the method of the invention 
comprises five steps referenced (a) to (e). 




namely ^i^(x,y) = 



XjVk - XkVj + (Yj - Vk )x + (Xk - Xj)y 
XjYk -XkYj +Xkyi -x^y, +x.y.^-x.^y, 



according to this method is given by the following linear system : 

D^^i = D^-[R^ + aJsNr^VE"^ c> - H.5D^"^ = VE^ (1) 

with : 

- D'^"'^ a column matrix of 2N elements including the components dx and dy 
5 of the nodal vectors on iteration k+1, N being the number of nodes of the 

meshing in the current step : 

- D'^ a column matrix of 2N elements including the components dx and dy of 
the nodal vectors on iteration k ; 

-H = [R^ + a.l2N] 
10 - I2N the identity matrix with the dimension 2N ; 

- VE^ = I column matrix of 2N elements in which 

N elements Vx.nE'* and N elements Vy.nE", 
n denoting a node of the meshing and taking in turn 
the values (1 ...N) ; with 

V„E>=2. X Y. DFD,(x,y) /'^"'"^''"'^'' .r(x,y) 

eesupp(n) (x,y)ee CX 



Vv,nE^=2- Z . Z DFDJx.y).^M^ll^^.^:(x,y) 



eesupp(n) (x,y)ee 



where DFDk (x,y) = Y2(x,y) - Yi (x-dx,y-dy) on iteration k where sup p(n) 
represents the support of the basic function T,^(x,y) attached to the node 
n, that is all the triangles having the node n for a vertex ; 



[^Rk,yx pjk.yy J 



-R^^l^,^!^^ r^kwl ^ square matrix with dimension 2N 

where 

RL^=2. Z £ P'-"'":^'^-''' Tx(x,y)X°(x,y) 

e£supp(tnn) (x,y)Ee CX J 



zf "'"'l"-'''f'''''l'^"'''' k(x,y).r(x,y) 

e.sut;(mn)(x7)^eL dX X dy j 

R..=2. E Z f "'"'l''"''' T'^:(x.yWx,y) 

eesupp(mn) (x,y)€e V_ Cy ^ 

where m and n denote nodes of the meshing and take in turn the values 
(1 ... N) and where supp(nm)=supp(n) o supp(m). 

- a = max,|v,E'|.||^„||) 

where p^„\\\s a functional norm of ^n- The two most advantageous norms 
are : 

!l'F„il = sup,.,>,.„„<„,|'r„(x,y)|= 1 or 



y |SUP p^njl (x,y)esupp(n) 

|supp(n)| denotes the cardinal number of the discrete region sup p(n). 

At the end of this optimisation phase, there are available N nodal vectors each 
relating to one node of the meshing. 

According to a variant embodiment of the adaptive gradient, It is possible to 
consider using a decomposition technique l<nown as an "LDL* profile" in technical 
language so as to resolve the linear system (1) and accelerate the treatment. This 
technique is described in the work entitled "Matrix numerical analysis applied to 
the art of the engineer" by Theodore Lascaux, Volume 1 , pp 295-299, published 
by Masson, 1986. 
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According to an advantageous characteristic of the invention, the meshing is next 
locally refined via the division of triangles when the variation between the image I2 
and the image I2 on these triangles is too high. This is why according to the step 
(c) of the method a variation E is calculated between the image I2 and the image 
5 I2 for each triangle e. The variation E is defined as follows : 

E= ;^DFD^(x,y) 

(x,y)ee 

With DFD(x,y) = Y2 (x,y) - Yi (x-dx.y-dy) 

Of course, so as to calculate this variation for each triangle, it is necessary to 
firstly have calculated the value of the movement vectors of all the points of the 
10 image I2 by means of interpolation from the nodal vectors obtained at the end of 
step (b). 

Then in accordance with step (d), the meshing is refined on a discrete fraction of 
all the triangles of the model. This fraction is determined according to a criterion 
15 relating to the variations E previously calculated in step (c). So as to carry out this 
refining, it is possible for example to classify the triangles of the model by a 
decreasing order of their variations E and subdividing the X first triangles of this 
classification into smaller triangles. X is a predetermined fraction of the number of 
finished elements in the model, for example half. 

20 

So as to locally refine the meshing, it is also possible to compare all the variations 
E calculated in step (c) with a threshold variation which depends on the size of the 
finished element in question and of subdividing into smaller finished elements the 
finished elements whose variations E are greater than the threshold variation. 

25 

The subdivision of a triangle e into four smaller triangles is shown on figure 2. The 
triangle e is defined by the three vertices Pi, P2 and P3 having for respective 
movement vectors d^, d^, dj . So as to subdivide it into four, three new nodes P4, 
P5, Ps are defined in the middle of the three sides P1P3, P1P2, P2P3 of the triangle. 
30 Allocated to each of these three new nodes is a movement vector equal to the 
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average of the movement vectors of the two tops of the side to which it belongs, 
respectively (H,+^3)/2, (^,+^2)^2, (d2 + 3*3)72 . The adjacent triangles to the 
triangle e whose side is P1P2, P2P3 or PiPsare then subdivided into two or three. 

5 Thus a model of finished elements is obtained whose meshing has been locally 
refined. According to step (e), steps (b), (c) and (d) are repeated at the end of the 
preceding step (d). This succession of steps is then repeated until a stoppage 
criterion is satisfied. This stoppage criterion is for example a predetermined 
number of finished elements to be reached at the end of step (d). 

10 

It is also possible to stop the method when the variations E of all these finished 
elements of the model obtained at the end of the preceding step (c) are lower 
than a threshold variation. 

15 According to an improved embodiment shown on figure 3, the steps (b) to (e) are 
carried out by depending on several resolution levels of images U and I2. The aim 
of this variant is to improve and accelerate the convergence of the calculations of 
the movement vectors. In order to achieve this, first of all for each pair of 
numerical images h and I2, a set of R images l[ is defined with a level of 

20 resolution r and luminance Y, , r taking in turn the values (R-1, R-2,...0) and i the 
values (1,2) and then steps (b) to (e) are carried out for each level of resolution r 
from the level of resolution r=R-1 to the level r=0. It is to be noted that the images 
1° and 1° correspond to the numerical images U and I2. 

25 In practice, the images Ij are obtained by the filtering of the image h using a linear 
low pass filter only allowing 1/2'' of the spectral band of the image in question in 
the directions x and y, that is a pulse response filter h[, having a pass-band 
BPr=[1/2''^\ 1/2^'"^) in the space of standardised frequencies [-1/2, 1/2]. The image 
jl" is defined by the following equation : 



30 



M M 

Y;(x,y)= 2 2:Y(x-u,y-v)hK 

u=-Mv=-M 
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The filter used is for example an approximation of an ideal filter and its pulse 
response is defined as follows : 



As indicated previously, this optimisation on several resolution levels is able to 
improve and accelerate the convergence of the calculations of the movement 
vectors. It is to be noted that the number of resolution levels R selected may differ 
from the number of successive refinings carried out on the meshing. 

According to a preferred embodiment, compactness constraints are added on 
each triangle of the model so as to prevent the triangles from reversing. 

The compactness of a triangle with vertices Pi, Pj and Pr is defined by the 
following equation : 



_r M 

with S- Ts' 



- M < n < M 



and s[ =2B.sinc(27iB^n) = 2B 



sin27tBrn 
27cB,n 



where B and M are natural integers 
M = +00 in the ideal case. 




with C(Pi,Pj,Pk) G ]0,1[ ; and S(Pi,Pj,Pk) et P(Pi,Pj,Pk) representing 
respectively the surface and perimeter of the triangle (Pj, Pj, Pk) 



If the compactness of a triangle is prevented from tending towards zero, it is also 
prevented from reversing. This is why, so as to avoid the reversals of triangles, 
each triangle must verify the following constraint : 
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C (Pi + , Pj + dp^ , P, + dp j > K X C {P, , Pj , P, ) 

K X C (Pi , Pj , P J - C (P; + dp, , Pj + dp^ , P, + dp J < 0 
og,(D)<0 ; e = triangle(Pi,Pj,Pj 

where K is a parameter fixing the authorised compactness variation and 
"D is the column vector of the movement vectors of the nodes of the model. 

According to the invention, this constraint is associated with each triangle at the 
time the movement vectors are optimised. The step for optimising the movement 
vectors amounts to a system of the type : 
mineE(ci) 

.g^(D)<0 Veel 
where: 

- E(^} represents the variation between the image b and the said image h ; 

- Qe is a constraint related to the triangle e ; 
- 1 is the set of the triangles of the meshing. 

So as to resolve the optimisation problems under constraints, the so-called 
increased Lagrangian technique is used. This technique is described in the work 
entitled "Theories and algorithms" by Michel Minoux, Volume 1, pp 257-260, 
published by Dunod 1983. This technique combines two optimisation techniques : 
Lagrangian optimisation and the optimisation of external penalties. 

According to this technique, resolving the preceding system amounts to resolving 
the system without constraints according to : 

min^(^E(D)+2:G(g,p)x„rjj 

Where re is a penality element. 
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^eis a Lagrange multiplier 

G is an increased Lagrangian detemnined by the equation : 

Kg,(D)+r,g,(^) si r,>0 et g,(^)>0 
G(g3(D)x„rJ=k,g,(^) si r,=0 et g,(^)>0 

[o si ge(D)<0 

The constraints ge have been previously linearised by the Taylor formula to the 
order 1 : 



(P,) 



The optimisation method is then the following : 

- k=0 is initialised 

- X = 0 and r=0 are placed, X g iR'^and r e SI"" , and m denotes the number 
of triangles of the mode! 

- then the minimum 5Dk+i(?^,r) is determined so that on iteration k+1, 

-H.5D'"^ =VE'-CV 
where Y-(?.,r), C* is a matrix of ^l^"" x9^"" 

CV forms a matrix of the linearised constraints having for coefficients the 
following values : 

0 si ge(^)<0 



sinon 



5a.^ge(Pi.Pj.Pj si P,=P,Pj ou P, 
0 sinon 
where e represents the triangle with the vertices Pj, Pj and Pk. 



Then X is updated by the Uzawa algorithm and r is increased. Then the 
preceding operation is repeated until all the constraints are verified before moving 
on to iteration k+2. It is to be noted that details of the Uzawa algorithm are given 
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in the work entitled "Theories and algorithms" by Michel Minoux, Volume 1, 
published by Dunod 1983. 

According to a final embodiment, it is also possible to introduce constraints so as 
to avoid the flowing over of the meshing obtained after applying the movement 
vectors beyond the range of the image li. This embodiment consists of forcing the 
peripheral nodes of the meshing to remain on the edges of the image after 
applying movement vectors. In order to do this, the abscissae components 50"*^ 
for the peripheral nodes on the left and right edges of the image b are cancelled 
on each iteration k. Similarly, the ordinate components 50"*^ for the upper and 
lower edges of the image b are cancelled on each iteration k. 

In the case of a coding application, associated with the meshing is a partially 
quaternary tree. On each subdivision of the meshing (step d), an additional level 
is added in the tree. Each level of the tree then represents a meshing level and 
each node of the tree represents a triangle of the corresponding meshing level. 
The binary train generated during the coding step is obtained via a reading of the 
tree level by level. In this binary train, the movement vectors associated with each 
node of the tree are advantageously coded differentially with respect to the 
movement vectors of their father node when the latter exists and are arranged 
level by level. The corresponding decoding step consists of regenerating this tree 
from the binary train received derived from the encoder. 
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CLAIMS 

1. Method for estimating the movement between two numerical images U and I2 
with luminance Yi and Y2 for generating for each point of coordinates x,y of the 

5 image b a movement vector d(x,y)=(dx,dy) so as to form an image I2 from the 
image I1 with luminance Y2(x,y)=Yi(x-dx,y-dy) which is an approximation of the 
image b, characterised in that it comprises the following steps : 

(a) defining an initial model of finished elements comprising a meshing 
whose nodes are points of the image b, a movement vector at each node of said 

10 meshing, and an interpolation formula for calculating the value of the movement 
vector of each point of the image I2 from the values of the movement vectors of 
the nodes of the mesh to which it belongs, 

(b) globally optimising the values of all the movement vectors of the model 
according to a differential method, 

15 (c) calculating a variation E between the image I2 and the image I2 for each 

finished element or mesh, 

(d) carrying out a finer meshing on a discrete fraction of all the finished 

elements determined according to a criterion relating to the variations E and 

allocating a movement vector to each new meshing node, 
20 (e) repeating the steps (b), ( c ) and (d) on the model of finished elements 

obtained at the end of the preceding step (d) until a stoppage criterion is satisfied. 

2. Method according to claim 1 , characterised in that, so as to carry out a finer 
meshing on a discrete fraction of all the finished elements in step (d), said set of 

25 finished elements is classified in the decreasing order of their variations E and the 
X first finished elements of this classification are subdivided into smaller finished 
elements, X representing a predetermined fraction of the number of finished 
elements of the set. 

30 3. Method according to claim 1 , characterised in that, so as to carry out a finer 
meshing on a discrete fraction of the set of finished elements in step (d), the set of 
variations E calculated in step (c) is compared with a threshold variation which 
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depends on the size of the finished element In question, and the finished 
elements whose variations E are greater than the threshold variation are 
subdivided Into smaller finished elements. 

4. Method according to one of claims 1 to 3, characterised in that said stoppage 
criterion is a predetermined number of finished elements constituting the model of 
finished elements to be reached at the end of step (d). 

5. Method according to one of claims 1 to 3, characterised in that said stoppage 
criterion of step (e) is satisfied when the variations E of the set of finished 
elements of the model obtained at the end of the preceding step (d) are smaller 
than a functional threshold variation which depends on the size of the finished 
elements in question. 

6. Method according to one of claims 1 to 5, characterised in that in addition for 
each numehcal image li and b, a set of R images ir with a level of resolution r 
and luminance Y; with r taking the values (0,....R-1) and I the values 1 and 2 Is 
defined, the images 1° and 1° corresponding to the numerical Images h and b, 
and in that the steps (b) to (e) are carried out for each resolution level r from the 
level r=R-1 to the level r=0. 

7. Method according to claim 6, characterised In that the sets of R images with 
resolution level r are obtained by filtering the Images li and h along the two 
directions x and y using a low-pass filter with a pulse response h^ , each image l[ 
being defined by the following equation : 

M M 

Yi (x,y)= 2 ZYi(x-u,y-vKh:, 

u=-Mv=-M 

with M a natural integer. 

8. Method according to claim 7, characterised In that the pulse response h^, is 
defined as follows : 
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gr M 

hn=— avec S=^s[, 



M 




s;, =2B.sinc(27cB,n) = 2B 



sin27iB,n 
2iiB,n 




B being a natural integer. 

9. Method according to one of the preceding claims, characterised in that the 
initial movement vectors are nil vectors. 

10. Method according to one of the preceding claims, characterised in that the 
variation E between the image b and the image b for each finished element e is 
defined by the following equation : 



where DFD(x,y) = Y2 (x,y) - Yi (x-dx.y-dy) 

1 1 . Method according to one of the preceding claims, characterised in that the 
interpolation formula for calculating the value of the movement vector of a point P 
of coordinates {x,y) in the image b belonging to the finished element e with 
vertices Pi, Pj and Pk with respective coordinates (Xi,yi), (Xj,yj) et (Xk,yi<) is the 
following : 

H(x,y)= 2^i'(x.y)-^(xi.y.) 

I=i,j.k 

where is a function of the form 
[T| (x,y) = a|+P|X + Yiy (x,y)ee 



E= J^DFD^i 



(x,y) 
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12. Method according to one of claims 1 to 11, characterised in that the 
differential method for optimising the movement vectors is the Gauss-Newton 
method. 

13. Method according to one of claims 1 to 11, characterised in that the 
differential method for optimising the movement vectors is the Marquardt 
extension of the Gauss-Newton method. 

14. Method according to one of the preceding claims, characterised in that a 
compactness constraint is imposed on each finished element at the time of 
optimising the movement vectors of the model of finished elements, said 
constraint consisting of preventing the compactness of each finished element 
from tending to zero. 

15. Method according to claim 14, characterised in that the compactness 
constraint on a finished element e with vertices Pj, Pj, Pk and compactness 
C {Pi,Pj,Pk) is defined by the following equation : 

C (Pi + 5. , Pj + , P, + dp J > K X C (Pi , Pj , P J 

where d^, et 5^ represent the movement vectors of the vertices Pi, 
Pj, Pk during the optimisation step, and K is a compactness parameter. 

16. Method according to claim 14 or 15, characterised in that the optimisation of 
the movement vectors under constraints on the finished elements is resolved by 
the increased Lagrangian technique. 

17. Method according to claim 16, characterised in that the constraints are used in 
a linearised form in the increased Lagrangian technique. 

18. Method according to one of claims 12, 13 or 17, characterised in that the 
methods for optimising the movement vectors use an LDL* profile technique. 
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19. Application of the method for estimating movement between two numerical 
images according to one of claims 1 to 18 for coding images, characterised in that 
the fractional subdivision of the meshing carried out in step d) of the movement 
estimation method is associated with a partially quaternary tree in which each 
level represents a meshing level and each node represents a triangle of the given 
level, and in that what is generated is a binary train describing said tree. 

20. Application of the method for estimating movement between two numerical 
images according to claim 19, characterised in that the movement vectors 
associated with each node of said tree are encoded differentially with respect to 
the movement vectors of their father node when the latter exists and are ordered 
in said binary train along a width passage of said tree. 

21. Application of the method for estimating movement between two numerical 
images according to one of claims 1 to 18 for decoding images, characterised in 
that the fractional subdivision of the meshing carried out in step d) of the 
movement estimation method is associated with a partially quaternary tree in 
which each level represents a meshing level and each node represents a triangle 
of the given level, and in that said tree is generated from a binary train of encoded 
data describing said tree. 

22. Application of the method for estimating movement between two numerical 
images according to claim 21, characterised in that the encoded data relating to a 
given level of the tree are collectively regrouped in the binary train so as to 
generate the tree level by level as the train is read. 

23. Application of the method according to one of the preceding claims to at least 
one of the ranges belonging to the group of the following ranges : 

- compression of sequences of images, and 

- compression of data in spaces larger than 2. 
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